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ABSTRACT 


The  flow  between  compartments  in  physical  and  biological  systems  is 
treated  as  a  special  case  of  a  more  general  theory  of  transitions  between 
any  two  distinct  sets  A, A.  Interest  is  focused  on  the  flow  rate  from 

each  set,  i.e.,  the  rate  at  which  elements  from  that  set  appear  in  the 
other;  and  on  the  entry  rate  from  each,  i.e.,  the  rate  at  which  elements 
from  the  set  leave  to  enter  the  region  not  part  of  either  set.  In 
particular,  the  two  flow  rates  are  completely  determined  by  means  of 
explicit  expressions  for  their  rates  (Theorem  l)  and  difference 
(Theorem  ll)  in  terms  of  the  two  entry  rates.  An  application  to  bio¬ 
logical  transport  problems  extends  a  result  of  Dantzig  and  Pace  [2]  by 
demonstrating  that  for  a  system  of  channels  each  narrow  enough  to  effect 
a  "lining-up"  of  particles,  counter-gradient  flows  may  result,  i.e., 
flows  for  which  the  flow  rate  is  greatest  from  the  compartment  with  the 
smallest  entry  rate. 


ON  STEADY-STATE  INTERCOMPARTMENTAL  FLOWS 


In  this  paper,  we  study  the  flow  between  compartments  in  physical  and 
biological  systems  as  a  special  case  of  a  more  general  theory  of  transitions 
between  any  two  distinct  sets  A,  A.  We  define  a  mesh  between  A,  A  to  be  any 
sequence  of  disjoint  sets 


A  “  S0’  Sl . Sn’  Sn+I  "  A 


with  the  property  that  an  element  of  any  S^  may  transition  to  other  sets  only 
by  a  sequence  of  forward  and  backward  movements  between  consecutive  sets. 

Elements  of  Sg  may  enter  the  i nter ior  of  the  mesh  by  transitioning  to  , 
and  the  rate  at  which  this  occurs  is  denoted  E  and  called  the  forward  entry 
rate.  Similarly,  elements  of  Sn+j  may  enter  the  interior  of  the  mesh  by 
transitioning  to  Sn>  and  this  backward  entry  rate  is  denoted  E.  Elements  which 
enter  the  interior  from  Sg  are  called  a-elements  until  they  leave  the 
interior  of  the  mesh  either  by  dropping  back  from  Sj  to  Sg,  which  occurs  with 
rate  R,  called  the  return  rate,  or  by  transitioning  from  Sn  to  Sn+j  which 
occurs  with  rate  F,  called  the  forward  flow  rate.  Between  entering  the  interior 
of  the  mesh  and  leaving  it,  either  by  dropping  beck  into  A  or  flowing  into  A, 
an  element  may  transition  by  a  sequence  of  forward  steps,  from  Sk  to  Sk+J ,  and 
backward  steps,  from  to  $k_j.  Similarly,  a-elements  may  enter  the  interior 

of  the  mesh  from  Sn+|  and  leave  either  by  dropping  back  from  Sn  to  Sn+j  , 
with  rate  R,  or  by  transitioning  from  S^  to  Sg,  with  rate  F. 


A  set  is  said  to  have  steady  content  if  its  content  of  a-elements  and  of 


•-elements  Is  independent  of  time. 

A  mesh  is  said  to  have  steady  content  if 

(I)  E-F  +  R.E-F  +  R 

•Iways;  i.e. ,  if  the  union  of  the  sets  has  steady  content. 

A  mesh  i s  cal  led  Markovian  i f  for  any  k"l , . . .  ,n,  the  probab  i  1  i ty  of  di rect ion 
of  movement  of  an  eJement  of  $k,  from  to  an  adjacent  set,  depends  only  on  the 

direction  from  which  it  entered,  and  not  on  any  other  aspect  of  its  past  history. 

For  such  meshes,  for  elements  which  entered  on  a  forward  step,  let  be 

the  probability  that  the  direction  of  its  next  movement  will  continue  to  be  forward, 
and  let  (l-ark)  he  the  probability  that  it  will  be  backward.  For  elements  which 
entered  on  a  backward  step,  let  3^  be  the  probability  that  the  direction  of 
its  next  movement  will  continue  to  be  backward,  and  let  (l-3k)  he  the  probability 
that  it  will  be  forward. 

A  Markovian  mesh  is  said  to  be  in  steadv-state  if  all  ^  ,  3^  are  independent 
of  time;  and  if  for  every  k  ■  l,...,n,  the  set  has  steady  content. 

We  can  now  prove: 

THEOREM  I:  For  Markovian  meshes  in  str>  '  -rtate.  the  ratio  of  flow  rates  between 
compartments  A  and  A  en  by  the  equation 


For  example,  if  A  and  A  represent  physical  compartments  connected  by  a  channel 
such  that  the  movement  of  elements  through  the  channel  is  unaffected  by  the 
presence  of  other  elements  in  the  channel,  then  for  any  mesh  that  might  be  defined 


by  subdividing  the  channel  it  follows  that  ■  0^  *  I  for  all  k;  and  hence  the 


flow  rates  are  proportional  to  the  entry  rates.  This  corresponds  to  Fick's  law 
for  free  diffusion.  On  the  other  hand,  as  in  [2]  and  [3],  the  channel  might 
instead  be  so  narrow  that  elements  line  up  in  a  single  file  of  n  elements,  and 
in  this  case  the  direction  of  next  movement  of  any  element  in  the  channel  will 
be  the  same  as  the  next  entry  into  the  channel  as  a  whole.  Thus  the  ratios 


°k  E  .  I-Bk 
'-“k  '  i  Bk  E 


a 

E  k  E 

—  ,  whence  "3“  ■  “  •  Hence  in  the  case  of  narrow  channels, 

7  K  7 


it  follows  from  the  theorem  that 


F 

F 


so  that  the  larger  n 


unidirectional  the  flow  becomes,  always  favorirg  the  direction  of  the 
rate. 


is,  the  more 
larger  entry 


The  proof  of  the  theorem  is  by  induction,  and  is  clear  for  n  ■  1.  We  assume 
the  theorem  true  for  n  ■  k-1,  and  apply  it  to  the  mesh  Sg  ,  Sj  ,...,S^  ,  defined 
by 

S'  ■  A  S'  ■  S  S'  ■  S  S'  ■  S  US  S'  *  A 

b0  A’  bl  V"*’V2  k-2  ’  k-1  Vl  *k  •  *k 

where  $k_j  U  Sk  denotes  the  union  of  the  two  sets  Sk_|  and  Sk*  This  mesh 
is  also  Markovian  and  in  steady-state.  We  obtain  by  the  inductive  assumption  that 


(2) 


F 

F 


a,...  a, 


ot 


ILL :  k-2  k-1 
V"  ek-2Pk-l 


f 


whe  re  a'  . 

k-1 


is  the  probability  that  an  element  which  transitions  from 


3 


to  will,  if  it  leaves  ■  S^j  U  S^,  next  transition  to  S^+j  ; 

and  is  the  probability  that  an  object  which  transitions  from  S^+j  to 

will,  if  it  leaves  S^_j  U  S^,  next  transition  to  S^*  It  is  easily  verified 
that  is  given  by  the  infinite  sum 

Vl^k  +  “k-l  ^’V^’^k-l^k  +-+Vl(,*ak)P(1'Pk-l)Pffk  + 
whence 


(3) 


Similarly,  we  obtain  that 


o1 

k-l 


oc  a 

k- 1  k 


,-(,-V(,'ek-.) 


(4) 


PhA 


k-l 


,"(1”ak)  (l”\-l  > 


and  the  theorem  is  immediate  upon  substitution  of  (3),  (4)  into  (2). 

In  physical  steady  content  situations,  we  are  often  interested  in  meshes 

that  have  the  t ransfer  property;  i.e.,  that  the  rate  E  at  which  a-elements 

enter  Sj  from  Sq  is  equal  to  the  total  rate  F+R  at  which  a-  and  a-elements 

leave  S  to  go  into  S  and  the  rate  E  at  which  a-elements  enter  S  from 

n  3  n+l  n 

Sn+|  is  equal  to  the  total  rate  F+R  at  which  a-  and  a-elements  leave  Sj  to 
go  into  Sq.  The  analytical  definition  of  this  property, 

E-F  +  R  ,  E-F  +  R 

taken  with  (l),  provides  an  immediate  characterization  of  these  meshes: 

THEOREM  II:  Any  mesh  with  steady  content  h?s  the  transfer  property  if  and  only  if 
(5)  F  -  F  -  E  -  E 

is  satisfied. 
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Since  any  mesh  in  steady-state  has  steady  content,  we  obtain  trivially  for  this 
important  case  the 

COROLLARY:  In  steadv-state .  a  mesh  satisfies  the  transfer  property  if  and  only 

if  the  difference  in  flow  rates  between  A  and  A  is  equal  to  the 
difference  in  entry  rates. 


Appl  icat ions: 

Some  investigators,  in  modeling  transport  in  biological  systems,  in 
particular  [2]  the  active  transport  of  sodium  and  potassium  ions  and  [3]  the 
exchange  of  potassium  t?nd  its  labelled  isotope  across  the  membranes  of  living 
mils,  have  postulaced  as  part  of  their  models  systems  of  channels,  each  channel 
ow  enough  to  effect  a  "lining  up"  of  particles.  Such  channel  s,  when  the  system 
is  physically  in  steady-state,  are  either  at  rest,  in  which  case  they  contain  a 
fixed  number  of  particles  ordered  in  some  way;  or,  they  are  in  a  transient  state 
caused  by  the  entry  of  a  particle  into  one  end  or  the  other  of  the  channel  from 
one  of  the  two  compartments  A,  A  between  which  transport  is  taking  place.  When 
an  entry  does  occur,  it  sets  off  a  chain  reaction  which  is  idealized  as  just  one 
particle  moving  at  any  time.  As  it  approaches  the  particle  which  at  rest  was  in 
the  position,  it  either,  with  probability  t^,  passes  it  and  proceeds  on  to 
the  next  particle;  or  it  st  r  i  kes  and  replaces  the  particle,  which  itself  then 
moves  on  to  the  next.  The  transfer  property  thus  holds  for  such  channels,  so  that 
the  difference  of  their  flow  rates  is  equal  to  the  difference  of  their  entry  rates 
(Corollary  to  Theorem  ll). 
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Fig.  I . 


Let  us  now  consider  the  ratio  of  the  flow  rates.  We  define  a  mesh  [See 
Fig.  1.]  between  A, A  by  conceptually  dividing  each  channel  into  n  consecutive 
sections  such  that  each  section  contains  exactly  one  particle  at  rest.  For  sim¬ 
plicity,  n  is  assumed  the  same  for  each  channel,  and  the  channels  are  assumed 
mutually  exclusive.  The  set  Sj  is  defined  as  the  union  of  first  sections  of  the 
channels,  i.e.,  those  adjacent  to  S^;  the  set  is  defined  as  the  union  of 

second  sections,  i.e.  those  adjacent  to  the  first  sections;  and  so  forth.  The 
mesh  is  Markovian;  and  in  fact  it  is  easy  to  see  that 
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<*k  -  tk  +  0-0  -h  .  B.  -  t  +  0-0  ~ 

K  K  K  E+E  *  E+E 


th 

because  a  particle  striking  and  replacing  the  k  particle  has  probabi  1  ity 


— 3  ,  when  it  is  next  struck,  of  moving  towards  the  (k+l)S  particle;  and 
E+E  ^ 

probability  — ”  of  moving  towards  the  (k-l)5  particle.  Hence  we  obtain  from 
E+E 

Theorem  I  that 


(6) 


£-£•  IT  ••= 

F  E  k-1 


E  +  t.  E 


E  +  t.  E 
k 


For  channels  so  narrow  that  the  particles  line  up  in  single  file  and  it  is  not 
possible  for  one  particle  to  pass  another,  it  follows  that  tk  *  0  for  all  k. 

In  this  case,  (6)  was  essentially  given  by  Hodgkin  and  Keynes  [3]  and  (5)  by 
Dantzig  and  Pace  [2 lousing  methods  quite  different  from  ours.  Bigelow  [l  ]  has 
obtained  a  generalization  of  their  methods  and  results  when  t^  is  the  seme 
non-zero  constant  for  all  k.  The  case  tk"l  corresponds  to  free  diffusion, 
as  noted  earlier.  The  principal  motive  for  such  an  approach  in  [2]  and  [l]  is  to 
demonstrate  that  flows  against  electrochemical  gradients  are  possible  within  the 
rather  simple  framework  of  these  models,  and  it  suffices  for  the  purpose  of  this 
paper  to  now  qualitatively  demonstrate  this  point  in  the  most  generrl  case. 


Counter  Gradient  Flows 

When  dealing  with  chemica-l  solutions,  it  is  necessary  to  consider  from  the 
same  compartment  entries  of  particles  of  different  types,  in  particular  different 
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chemical  species.  For  each  species  i  we  then  have  "entr  E.  from  A 

and  E.  from  A,  with  the  total  entry  rates  being  given  by 

E  •  E  E.,E  -  L  E.  . 

i  ;  i 

i 

Further,  we  may  also  define  "flow  rates"  F.,  F.  for  each  species,  such  that 

F  -  I  F.  ,  F  -  I  F.  . 

i  .  i 

i  i 

If  we  assume  that  within  the  mesh  the  behavior  of  an  element  is  independent  of  its 

species,  then  the  flow  of  a  particular  species  F^  is  in  the  same  ratio  to  total 

flow  as  E  is  to  total  entries.  Thus 
s 

F  E  F  E 

e  »  »  »  • 

_*  ■  — * —  ’  — *“  ■  '■*■■■  * 

IF.  IE.  " 

i  i  IF. 

i 

and  this,  taken  with  (6),  yields  at  once  the  ratio  of  flow  rates 


(7) 


F 

F 

s 


n 


TT 

k-1 


Hlki 


a 


As  noted,  in  this  development  t^,  the  probability  of  passing  the  particle  in 

k*^  position,  is  assumed  independent  of  species.  A  more  general  model  is  under 

development  where  t^  depends  on  species. 

We  say  a  counter-Qradient  flow  exists  for  a  species  i  if  the  sign  of 

(F.-F.)  is  opposite  to  that  of  ( E .  —  E . ) .  In  other  words,  if  E.  >  E.,  a  counter- 
i  i  i  i  ii 

gradient  flow  will  be  said  to  exist  if  F.  <  F.;  and  analogously  for  the  case 
E.  <  E.. 

i  i 

Let  us  consider  a  particular  species  i«s  for  which  Es  >  E^  and  investigate 
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the  conditions  for  a  counter-gradient  flow  to  exist  for  this  species.  We  first 
note  thtft  if  the  total  flow  E  >  E,  then 

E  +  tkE  ^  E  +  tkE  ^ 

follows  for  0  £  t.£  I  ;  and  thus  it  follows  from  (7)  that  >  ' *  so  that  there 

k  F 

s 

is  no  counter-gradient  flow  for  the  species  s.  To  paraphrase,  we  have  shown  that 
for  the  case  E&  >  E^  ,  a  counter-gradient  flow  can  possibly  exist  only  when  the 
total  entries  by  all  species  E  >  E. 


We  thus  address  ourselves  to  the  question,  how  must  we  adju  st  the  rat io  of 

total  entry  rates  r  ■  ~  to  induce  a  counter-gradient  flow  for  the  species  s, 

E 

where  E^  >  E^  ?  By  definition,  a  counter-gradient  flow  exists  for  s  if  and 

only  if  <  I ;  or  noting  (7)  if  and  only  if 
F 

s  _ 

E 

(8)  i  >  t*  >  V**V 

s 


where  we  have  set 


E  +  fckE  ■  r  *  *k 

E  +  t  E  1  +  rtk 
k 


We  assume  that  the  probability  tk»  O^t^l  ,  of  bypassing  a  prrticle  in  the 

passage  is  independent  of  the  flow  rates  E,E.  It  is  easy  to  veiify  that  \k  is 

continuous  and  monotonical I y  non-decreasing  in  r  ■  ~  ,  and  ranges  between  its 

E 

extreme  values  I  and  tk  as  r  ranges  between  1  and  0.  Thus  it  follows  at 
once  that  the  product  is  continuous  and  monoton i cal  1 y  non-decreasing 

in  mr;  and  rrnges  between  its  extremes  1  and  t  ^ .  t n  as  r  ranges  between  1 
and  0,  whence 
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c 

THEOREM  III:  A  necessary  and  sufficient  condition  that  the  ratio  r  ■  ~  can  be 

r- 

C 


adjusted  so  that  a  counter-gradient  flow  exists  for  the  species  s. 
E 

where  ^  <  1 ,  is  that 
s 


(9) 


E 

>  7s-  >  t, . .  .t  . 

E  I  n 

s 


E 

(Note  that  the  case  >  1  >5  easily  handled  with  this  theorem  by  an  elementary 

s 

notational  exchange  of  barred  and  unbarred  quantities).  The  necessity  of  the  con¬ 
dition  follows  from  (8)  and  the  relation  just  obtained,  that  is  X . . . X  ^  t,...t  . 

Inin 

To  show  sufficiency,  suppose  (9)  holds.  Because  X(...Xn  ranges  continuously 

from  1  to  tj...t^f  there  always  exist  by  the  intermediate  value  theorem, 

X.  .... ,X  such  that  (8)  holds. 

I  n 

For  the  extreme  case  of  narrow  channels,  for  example,  t^*0  for  k,  so 

that  (9)  always  holds.  In  that  case,  it  is  easy  to  see  that  counter-gradient 

E  n  E- 

flows  will  exist  if  r  ■  “  is  chosen  so  that  r  <  T2-.  At  the  other  extreme, 

E  s 

for  "open"  channels  (i.e.  free  diffusion)  where  all  t^  ■  I,  condition  (8)  never 
holds  and  hence  counter-gradient  flows  can  never  occur. 
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